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Abstract 

We study the Korteweg-de Vries-type equation dtu = — d x (d x u + f(u ) — b(t, x)u ), where b is a small and bounded, 
slowly varying function and / is a nonlinearity. Many variable coefficient KdV-type equations can be rescaled into 
this equation. We study the long time behaviour of solutions with initial conditions close to a stable, 6 = 0 solitary 
wave. We prove that for long time intervals, such solutions have the form of the solitary wave, whose centre and scale 
evolve according to a certain dynamical law involving the function b(t,x), plus an Pf 1 (R)-small fluctuation. 


1 Introduction 

We study the long time behaviour of solutions to a class Korteweg-de Vries-type equations, with an 
additional term b(t, x)u. These equations, from now on called the bKdV, are of the form 

d t u = -d x (dlu + f(u) - b(t, x)u) , (1) 

where 6(f, x) is a real valued function and / is a nonlinearity. In this paper we consider a restricted class of 
nonlinearities. In particular, for monomial nonlinearities, we give a result only for f(u) = vr, corresponding 
to the modified KdV (mKdV). When 6 = 0, Equation Q reduces to the generalized Korteweg-de Vries 
equation (GKdV) 

d t u = -d x {dlu + f(u)). (2) 

A remarkable property of the GKdV is the existence of spatially localized solitary (or travelling) waves, 
i.e. solutions of the form u = Q c {x — a — ct ), where a £ R and c in some interval I. When f(u ) = u p and 
p > 2, solitary waves are explicitly computed to be 

Q c (x) = cp^Q{c^x), 
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where 



It is generally believed that an arbitrary, say H 1 ( R), solution to equation 0 eventually breaks up into a 
collection of solitary waves and radiation. A discussion of this phenomenon for the generalized KdV appears 
in Bona [Si. For the general, but integrable case see Deift and Zhou m- 

The mKdV equation is fundamental in many areas of applied mathematics ranging from traffic flow to 
plasma physics (see |23 El HUES) and arises from an approximation of a more complicated systems. The 
effects of higher order processes can often be collected into a term of the form b(t, x)u. Our main result 
stated at the end of the next section gives, for long time, an explicit, leading order description of a solution 
to the bKdV initially close to a solitary wave solution of the GKdV. 

We assume that the coefficient b and nonlinearity / are such that m has global solutions for H 1 ( K.) 
data and that m with b = 0 possesses solitary wave solutions. Precise conditions will be formulated in the 
next section. Here we mention that the literature regarding well-posedness of the KdV (b = 0, f(u) = u 2 ) 
is extensive and well developed. The Miura transform (see EH) then gives well-posedness results for the 
mKdV. Bona and Smith jSJ proved global wellposedness of the KdV in H 2 ( R). See also m, Kenig, Ponce, 
and Vega m have proved local wellposedness in H s ( K.) for s > — | and similar results are available for 
the generalized KdV ( b = 0, monomial nonlinearity f(u) = u p with p = 2, 3,4) [231. In particular, local 
well-posedness for the mKdV in H S (M .) with s > j and global well-posedness for s > 1 are known. More 
recently, results extending local wellposedness in negative index Sobolev spaces to global wellposedness have 
been proven EH0. There is little literature on global well-posedness of the bKdV in energy space, however, 
under a smallness assumption on the coefficient 6, Dejak and Sigal m proved global well-posedness in H 1 (R) 
of the bKdV with f(u) = u p , p = 2,3,4. They used results of m, and perturbation and energy arguments. 

Soliton solutions of the KdV equation are known to be orbitally stable. Although the linearized analysis 
of Jeffrey and Kakutani m suggested orbital stability, the first nonlinear stability result was given by 
Benjamin [li. He assumed smooth solutions and used Lyapunov stability and spectral theory to prove his 
results. Bona 0 later corrected and improved Benjamin’s result to solutions in ff 2 (R). Weinstein m used 
variation methods, avoiding the use of an explicit spectral respresentaion, and extended the orbital stability 
result to the GKdV. More recently, Grillakis, Satah, and Strauss m extended the Lyapunov method to 
abstract Hamiltonian systems with symmetry. Numerical simulations of soliton dynamics for the KdV were 
performed by Bona et al. See HElllIZI. 

For nonlinear Schrodinger and Hartree equations, long-time dynamics of solitary waves were studied 
by Bronski and Gerrard nni, Frohlich, Tsai and Yau El, Keraani m, and Frohlich, Gustafson, Jonsson, 
and Sigal pi E2. For related results and techniques for the nonlinear Schrodinger equations see also 

im 03123133133 HD 03133133. 

In our approach we use the fact that the bKdV is a (non-autonomous, if b depends on time) Hamiltonian 
system. As in the case of the nonlinear Schrodinger equation (see El), we construct a Hamiltonian reduction 
of this original, infinite dimensional dynamical system to a two dimensional dynamical system on a manifold 
of soliton configurations. The analysis of the general bKdV immediately runs into the problem that the 
natural symplectic form to is not defined on the tangent space of the soliton manifold. In this paper we prove 
the main theorem in the cases where the symplectic form is well defined on the tangent space. One such case 
is when the nonlinearity is f(u) = u 3 . For the general case see H2. We remark here that the dynamics for 
the special case considered here include the higher order correction terms for the scaling parameter c, which 
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cannot be included in the general case. 


Acknowledgements 

We are grateful to I.M. Sigal for useful discussions. 


2 Preliminaries, Assumptions, Main Results 


The bKdV can be written in Hamiltonian form as 


d t u = d x H' b [u), (3) 

where H' b is the L 2 (R) function corresponding to the Frechet derivative dH b in the L 2 (R) pairing. Here the 
Hamiltonian H b is 


r°° 1 1 

H b (u) := / ~{d x u) 2 — F(u) + -b(t, x)u 2 dx, 

J — OO " " 

where the function F is the antiderivative of / with F(0) = 0. The operator d x is the anti-self-adjoint 
operator (symplectic operator) generating the Poisson bracket 

1 r°° 

{G u G 2 } := - / G\ {u)d x G', 2 (u) - G' 2 (u)d x G[(u) dx, 
defined for any Gi, G 2 such that G\, G' 2 € ffl(R). The corresponding symplectic form is 

u{v i, v 2 ) ■■= - / vi{x)d~ 1 v 2 {x) - v 2 {x)d x v\(x) dx, 

defined for any vi,v 2 £ L 1 (R). Here the operator d~ l is defined as 

d~ l v{x) := / v{y)dy. 

J —OO 

Note that d■ d x = I and, on the space {u £ L 2 (R) | f^° udx = 0}, d~ x is formally anti-self-adjoint with 
inverse d x . Hence, if Jv\(x) dx = 0, then ui{v i, v 2 ) = vi(x)d~ 1 v 2 (x ) dx. 

Note that if b depends on time t, then equation @ is non-autonomous. It is, however, in the form of a 
conservation law, and hence the integral of the solution u is conserved provided u and its derivatives decay 
to zero at infinity: 

d r t 
~ / udx = 0. 
dtj_ oo 

There are also conserved quantities associated to symmetries of JD when 6 = 0. The simplest such corre¬ 
sponds to time translation invariance and is the Hamiltonian itself. This is also true if b is non-zero but time 
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independent. If the potential 6 = 0, then 0 is also spatially translation invariant. Noether’s theorem then 
implies that the flow preserves the momentum 

P(u) := |ll^lli- 

In general, when 6 ^ 0 the temporal and spatial translation symmetries are broken, and hence, the Hamil¬ 
tonian and momentum are no longer conserved. Instead, one has the relations 

d 1 r°° 

-H b (u) = -J ( d t b)u 2 dx , (4) 

s f( " ) = 2 LJ' ud!C ' <5) 


where b'(t,x) := d x b(t,x). For later use, we also state the relation 

d 1 / ’°° 
dt 2 


bu 2 dx = f ^ u 2 d t b + b' f uf(u) — d x u ) 2 — F(u) j — b"ud x udx. 

) J — OO ^ \ " / 


( 6 ) 


Assuming 0 is well-posed in if 2 (R), the above equalities are obtained after multiple integration by parts. 
Then, by density of H 2 { R.) in U 1 (R), the equalities continue to hold for solutions in H 1 ( R). To avoid these 
technical details, we assume the Hamiltonian flow on H 1 ( R.) enjoys 0, 0 and @. 

Consider the GKdV, i.e. equation ©. Under certain conditions on /, this equation has travelling wave 
solutions of the form Q c (x — ct), where Q c a positive H 2 (WL) function. Substituting u = Q c (x — ct ) into the 
GKdV gives the scalar held equation 

-&lQ c + cQ c - f(Q c ) = 0. (7) 


Existence of solutions to this equation has been studied by numerous authors. See El 0. In particular, 
Berestyki and Lions El gi ye sufficient and necessary conditions for a positive and smooth solution Q c to 
exist. We assume g := — cu+ f(u) satisfies the following conditions: 


1. g is locally Lipscliitz and g(0) = 0, 

2. x* := inf{s > 0 | j* g{y) dy} exists with x* > 0 and g{x*) > 0, and 

3. lim s ^o p < — m < 0. 


Then, as shown by Berestycki and Lions, 0 has a unique (modulo translations) solution Q c G C 2 for 
c in some interval, which is positive, even (when centred at the origin), and with Q Cl d x Q c , and d 2 Q c 
exponentially decaying to zero at infinity ( d x Q c < 0 for x > 0). Furthermore, if / is C 2 , then the implicit 
function theorem implies that Q c is C 2 with respect to the parameter c on some interval Iq C M+. We 
assume that x m d™Q C £ L 1 (R.) for n = 1,2, m = 0,1, 2 so that integrals containing d™Q c are continuous and 
differentiable with respect to c. We also make the assumption that 



d c Q c dx = 0 


for all c £ I. This implies that 


/ ;X nx 

d c Qc(z) dz, / d 2 c Q c (z ) dz £ Lr 

-oo J — OO 


( 8 ) 

( 9 ) 
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To see this use the isometry property of the Fourier transform and the decay properties of d c Q c . The above 
requirements of Q c are implicit assumptions on the nonlinearity / and are true when f(u) = u 3 . Assumption 
© is a very important and restrictive requirement; it does not hold when f(x) = x p and p / 3. For the 
case where m does not hold see m 

The solitary waves Q c are orbitally stable if <5'(c) > 0, where <5(c) = P(Q C )- See Weinstein m the first 
proof for general nonlinearities. Moreover, in eh, Grillakis, Shatah and Strauss proved that S'(c) > 0 is a 
necessary and sufficient condition for Q c to be orbitally stable. In this paper, we assume that Q c is stable 
for all c in some compact interval I C Jo, or equivalently that S'(c) > 0 on I. For f(u) = u p , we have 
S'(c) = ||Qc=i||j, 2 ) which implies the well known stability criterion p < 5 corresponding to subcritical 

power nonlinearities. 

The scalar field equation 0 for the solitary wave can be viewed as an Euler-Lagrange equation for the 
extremals of the Hamiltonian H b = o subject to constant momentum P(u ). Moreover, Q c is a stable solitary 
wave if and only if it is a minimizer of H b —o subject to constant momentum P. Thus, if c is the Lagrange 
multiplier associated to the momentum constraint, then Q c is an extremal of 

A ca (u) := H b = 0 (u) + cP(u) (10) 

r°° i i 

= I ~(d x u ) 2 + -cu 2 - F(u)dx, 


and hence A' ca (Q c ) = 0. 

The functional A ca is translationally invariant. Therefore, Q ca (x) := Q c (x — a) is also an extremal of A co , 
and Q c {x — ct — a) is a solitary wave solution of 0 with 6 = 0. All such solutions form the two dimensional 
C°° manifold of solitary waves 

M s := {Qca | c £ /, a G R}, 
with tangent space Tq cil M s spanned by the vectors 

Cco : = daQca = -dxQca and : = d c Qca, (11) 

which we call the translation and normalization vectors. Notice that the two tangent vectors are orthogonal 
in L 2 (R). 

In addition to the requirement on 6 that 0 is globally wellposed, we assume the potential b is bounded, 
twice differentiable, and small in the sense that 

I W&| < eae?e?, (12) 

for n = 0,1, m = 0,1, 2, and n + m < 2. The positive constants e a , e x , and et are amplitude, length, and 
time scales of the function 6. We assume all are less than or equal to one. 

Lastly, we make some explicit assumptions on the local nonlinearity /. We require the nonlinearity to 
be k times continuously differentiable, with /^ bounded for some k > 3 and /(0) = /'(0) = 0. These 
assumptions ensure the Hamiltonian is finite on the space 7L 1 (R) and, since Q c decays exponentially (see 
0), both f(Q c ) and f'(Q c ) have exponential decay. 

We are ready to state our main result. Recall that Iq C R+ is an interval where Q c is twice continuously 
differentiable. 
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Theorem 1. Let the above assumptions hold and assume S'(c ) > 0 for all c in a compact set I C Iq. Assume 
e a < 1- Then, if e x < 1, eo and e t are small enough, there is a positive constant C such that the solution to 
m with an initial condition uq satisfying infQ^gMs ||wo — QcciWh 1 < e o can be written as 

u(x,t) = Q c (t){x - aft)) +£{x,t), 

where ||£(f)||#i = O ^eo + (e a e x eo)5 +c x +eA for all times t < Moreover, during this time 

interval the parameters aft) and c(t) satisfy the equations 

^ +6 '( a )^^ ( l) + 0 (( e o + ex + e «) 2 + (eaexeo)^(ex + et + eo)) , 

where c is assumed to lie in the compact set I. 

Sketch of Proof and Paper Organization. To realize the Hamiltonian reduction we decompose functions in a 
neighbourhood of the soliton manifold M s as 

'a — Q C a T C 

with £ symplectically orthogonal to Tq cii M s , i.e. £J -d~ 1 TQ ca M s . We show that there is an eo > 0 such that 
if the solution u satisifes the estimate infg ca ||u — Q C a\\m < to, then there are unique C 1 functions a(u) and 
c(u) such that u = Q c ( u )a(u) + £ with t;Ld~ 1 TQ ca M s . 

With the knowledge that the symplectic decomposition exists, we substitute u = Q ca + 1; into the bKdV 
0 and split the resulting equation according to the decomposition 

L 2 (R) = d~ l T Qca M s © 

to obtain equations for the parameters c and a, and an equation for the (infinite dimensional) fluctuation 
£. In Section 2] we isolate the leading order terms in the equations for a and c and estimate the remainder, 
including all terms containing f. In Sections [S] and 0 we establish spectral properties and a lower bound of 
the Hessian A” a on the space ( df 1 TQ ca M s ) . 

The proof that ||Cliff 1 is sufficiently small is the final ingredient in the proof of the main theorem. The 
remaining sections concentrate on proving this crucial result. We employ a Lyapunov method and in Section 
0we construct the Lyapunov function r c and prove an estimate on its time derivative. This estimate is later 
time maximized over an interval [0, T], and integrated to obtain an upper bound on T c involving the time T 
and the norms of f. We combine this upper bound with the lower bound on r c following from the results of 
Section Q and obtain an inequality involving HCIIff 1 - In Section [21 we solve the inequality to find an upper 
bound on HCHff 1 provided (IC(O)llff 1 is small enough. We substitute this bound into the bound appearing in 
the dynamical equation for a and c, and take e a e x and eo small enough so that all intermediate results hold 
to complete the proof. □ 


3 Modulation of Solutions 

As stated in the previous section, we begin the proof by decomposing the solution of 0 into a modulated 
solitary wave and a fluctuation £: 

u{x,t) = Q c (t)a(t){x) + f(x,t), (13) 
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with a, c, and £ fixed by the orthogonality condition 



( 14 ) 


where 


d x 1 : 9 i-> / ff(z) dz 


J —oo 


Note that d x 1 TQ ca M s is a subset of L 2 ( R) (see 0). 

The existence and uniqueness of parameters a and c such that £ = u — Q ca satisfies follows from 
the next lemma concerning a restriction of df 1 and the implicit function theorem. 

The restriction K of d~ l to the tangent space Tq ca M s is defined by the equation KP T = Prd x 1 Pt, 
where Pt is the orthogonal projection onto Tq c!l M s . In the natural basis {ClaiQa} of the tangent space 
Tq,_„ M s . the matrix representation of K is IV^ 1 fl ca , where 



N 


and 



(15) 


Recall that S(c) = |||Qc|| 2 2 . 

Lemma 1. If S' (c ) > 0 on the compact set I C R+ ; then the matrix fl ca is invertible for all c £ I, and 



(16) 


Clearly, || < |y^J 1 , where := infj S'(c). 


Proof. The lemma follows from the relations «= 0, = 0 and 

(Qa,Qc) = S'(c). 


□ 


Given e > 0, define the tubular neighbourhood U e := {m £ L 2 (R) | inf( c> a )e/ X R ||« — Qca||i 2 < e} of the 
solitary wave manifold M s in L 2 (R). 


Proposition 1. Let I C R+ be a compact interval such that c i—> Q ca is C' 1 (J). Then there exists a positive 
number £o = Eq(I ) = O (L<5'_| 2 ) dependent on I and unique C 1 functions a : U eo —* R_|_ and c : U Eo —> I, such 
that 



for all u £ U Eo . Moreover, there is a positive real number C = C{I) such that 


11^ Qc{u)a{u) II// 1 — f-' ^ inf M ||^ Qca II/7 1 


(17) 


for all u £ U Eo H H 1 ( R). 
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Proof. Let g := {g 1 , g 2 ) T £ R + x I and define G : R+ x / x H 1 ( R.) —> K 2 as 


G : (g, u) 


(Q ca -u,n ca C a ) \ 
(Q ca -u,n ca o )' 


where a = g 1 and c = g 2 . The proposition is equivalent to solving G(g(u),u) = 0 for a C 1 function 
g. Let go = ( ac) T . If G is C 1 , G(go,Q ca )=0, and d^F{g o,Qca) is invertible, then the implicit function 
theorem asserts the existence of an open ball B eo (Q ca ) of radius £q with centre Q ca , and a unique function 
9Qca ’■ Bs(Q ca ) R+ x 1 , such that G(gQ ca (u),u) = 0 for all u £ B eo (Q ca ). The first two conditions are 
trivial, and the third follows from Lemma |T| since d^ l G(go 1 Qca) = O ca . The radius of the balls B E (Q ca ) 
depend on the parameters c and a. To obtain an estimate of the radius, and to show that we can take £ 
independent of the parameters c and a, we give a proof of the existence of the above function gQ ca for our 
special case using the contraction mapping principle. 

We wish to solve G(g, u) = 0 for g := (g 1 , g 2 ) T with u close to Q ca in L 2 (R). Expand G(g, u) in g about 
= (ac) T : G(g,u) = G(g 0 ,u) + d^G{go,u)( y g - g 0 ) + R(l^,u), with R{g,u) = O (||/z - g 0 \\ 2 ) (G is G 2 ). 
Thus, we must solve g = /io — [9 M G(/uo,u)] _1 (G(no,u) + R(g,u)) for g. Clearly, since (9 M G(/uo, u) = f l ca , /i 
must be a fixed point of 

Hu^o (m) • Mo [G(mo> + R(/J>i w)]’ 

We now show that H ufl0 is a strict contraction, and hence has a fixed point. By the mean value theorem 


l|iWM 2 ) - ^o(Mi)ll < sup \\dvHuto\W\H 2 - Mill, 


where the supremum is taken over all allowed parameter values. Furthermore, we have 
d^H UIJ , 0 (g) = [d^G(g, u) - d u G(no,u)] 

= -Slfa [d»G{g, u) - d^G^g, Q C a) + d, l G{g, Q ca ) - d^G(g 0 , Q ca ) + 9^G(/x 0 , Q ca ) - d u G(g 0 , u)] 


Using the mean value theorem again, we compute that 


II^GM -a M G(/xo,«)|| < Gk5 + G 2 £ 

for some constants Gi and G 2 if ||/z — MoII < <5 and \\ u ~ Qcq||l 2 < £ - Combining all the estimates gives 

||ff«Mo(M 2 ) - Hullo (Mi)II < sup ||U“ a 1 || (Gk 5 + G 2 £) 11M2 Mi II • 

Thus, if S = \(C\ sup HO" 1 H)- 1 and e = fiC 2 sup ||fl^ a 1 ||)- 1 , then H u ^ 0 is a contraction. 

We now choose S and e so that maps Bg(g 0 ) to Bs(g o)- We have that 

11^0 - Moll < HO” 1 (G(no,u) + R(g,u)) II < supllO" 1 !! (|| G(hq,u) - G{g 0 ,Q C a) || + 0 (6 2 )). 

By the mean value theorem |j G(g 0 ,u) — G(g 0l Q co )|| < G 3 £. Thus, if we take <5 = 0 (sup ||0^ a 1 ||- 1 ) so that 
O (<5 2 ) < \ (sup HO- 1 !!) -1 S, then 


ll^o - Moll < G 3 sup ||O ca ||e+ -6. 


We now take e < \ (C 3 sup ||f2 ca 1 ||) 1 S to obtain || H u ^ 0 — /r 0 || < To complete the argument, take S to be 

the smaller of \ (Ci sup ||f 2 “ Q 1 ||) 1 and the above choice, and then e to be the smaller of j(C 2 sup ||^ 7 a 1 ||) _1 
and 8{AC 3 sup ||^^ a 1 ||) _1 - Using the bound on ||f2“ a 1 || we find that 

£ = 0(|/J 2 ) 

if sup ||f 2 ^ 0 1 || > 1 , or equivalently, when is sufficiently small. 

The above argument shows that there exists balls {B e (Q ca ) \ a G R+,c G 1} with radius e dependent 
only on the compact set I. Then, defining U Eo = U{B Eo (Q ca ) \ a G R+, c G 1} and pasting the C 1 functions 
9Qca together, into a C 1 function <?/ : U Eo —> R+ x I, proves existence of the required C 1 functions a(u ) and 
c(u). Uniqueness follows from the uniqueness of the functions gQ ca . 

Let u G U E , c G /, and a G R, and consider the equation 

^ Qc(u)a(u) ^ Qca T Qca Qc(u)a{u)‘ 

Clearly, inequality iTHT) will follow if || Q C a — Q c (u)a(u) | h 1 < C||m — QcciWh 1 f° r some positive constant C. Since 
the derivatives d c Q ca and d a Q C a are uniformly bounded in iJ 1 (R) over I x R, the mean value theorem gives 
that ||Qca — Qc(u)a(u) II h 1 < U||(c,a) T — (c(rt), a(u)) T ||, where the constant C does not depend on c, a. The 
relations gi(Q C a) = (c,a) T and gi(u) = ( c(u),a(u)) T then imply ||Q cq - Q c {v.)a(u) Wh 1 < C\\gi{Q ca ) - gi{u)\\. 
Again, we appeal to the mean value theorem and use the properties of U ca and that d u gi = d Pl G~ l d u G is 
uniformly bounded in the parameters c and a to obtain 11171) . □ 


4 Evolution Equations for Parameters a and c 


In Section 0 we proved that if u remains close enough to the solitary wave manifold M s , then we can write 
a solution u to m uniquely as a sum of a modulated solitary wave Q ca and a fluctuation £ satisfying the 
orthogonality condition O- Thus, as u evolves according to the initial value problem CD, the parameters 
a(t) and c(t) trace out a path in R 2 . The goal of this section is to derive the dynamical equations for the 
parameters a and c, and the fluctuation £. We obtain such equations by substituting the decomposition 
u = Qca + £ into JU and then projecting the resulting equation onto appropriate directions, with the intent 
of using the orthogonality condition on £. 

From now on, u is the solution of QJ with initial condition uq satisfying eo := inf Q ca eM s ||rto — QcuWh 1 < 
So, and To = Tq(uq) is the maximal time such that u(t) G U E for 0 < t < To. Then for 0 < t < To, u can be 
decomposed as in m and Cil- 

Proposition 2. Assume 8'(c) ^ 0. Say u = Q ca + £ is a solution to 171) . where £ satisfies \14\l - Then, if 
II^Hjji is small enough, e x < 1, and c G I, 


a 

c 


c — b(t, a) 

0 


+ b'(t, a) 


6(c) 

5’(c) 


0 

1 


+ Z(a, c, £), 


(18) 


where Z(a, c, £) = O (e a e 2 + e a ex\\S,\\H^ + llsll^i)- 
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Proof. Recall that the solitary wave Q ca is an extremal of the functional A ca . To use this fact we rearrange 
definition 113 of A ca to write the Hamiltonian H b as 

1 r°° 

H b (u) = A ca {u) - cP(u) + - / bu 2 {x) dx, 

^ J — OO 

where for notational simplicity we have suppressed the space and time dependency of b. Substituting Q ca + £ 
for u in m and using the above expression for H b gives the equation 

®Cca + cC" a + £ = 9 x A' ca (Q ca + £) — cd x [Q ca + £] + d x [{Q ca + £)&], 

where dots indicate time differentiation. Let Cq := A" a (Q ca ), 

5b := b(t, x ) — b(t, a) 


and 


5 2 b := b(t, x) — b(t, a) — b'(t , a){x — a). 

Taylor expanding A' ca (Q ca + £) to linear order in £, using that Q ca is an extremal of A ca and the relation 
Ca = -9 x Q ca gives that 

t = d x [{C Q + 5b + b{a) - c)£] + d x N' ca {£) - [a - c + b{a)]£ r a - cQ a 

+ b'{a)d x [{x - a)Q ca ] + d x [5 2 bQ ca \. (19) 

The nonlinear terms have been collected into iV' a (£) given by 123 in Appendix 1X1 

Define the vectors Ci : = Cm and C 2 : = Cm• Projecting U3 onto d x 1 Ci for i = 1,2 and using the 
antisymmetry of d x gives the two equations 

[d — c + 6(a)] [{C r a,d- 1 Ci) + (£,C»)] +c(Cm>'9x 1 Ci) + (C.&^C*) -d(£,Ci) = - b'{t,a)({x - a)Q ca , Ci) 

- (6 2 bQ ca , Ci) - (<»£, Ci) - (^Vca(C), Ci) - <£q£, Ci). (20) 

We can replace the term containing £ since the time derivative of the orthogonality condition {£,d x Q) = 0 
implies (£,9“ 1 Ci) = d(£,C i) — c(£,9 c 9“ 1 Ci)- Note that we have used the relation d a (i = —d x (i. Thus, in 
matrix form, nu becomes 


where 


{i + B)n c 


a — c+ b{t, a) 
c 


= X + Y , 


X := 


Y := - 


{6b^Cc r a) + {Ka(0,Cc r a) + 

(6%, Cm) + (KaiOXZa) + <£ Q £, Qa) 


( 21 ) 
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and 


B := 


(CO <£,C> A 0-1 
<£,C”> -(^dcd-^a) J ca 


We have explicitly computed ((x — a)Q ca , Q) to obtain the above expression for X. 

We now estimate the error terms and solve for a and c. The assumption on the potential implies the 
bounds 


|<5&| < £a£x(x — a) and |<5 2 6| < e a e 2 (:r — a) 2 . (22) 

Thus, Holder’s inequality and exponential decay of Q ca imply 

X = -b'(t,a)S(c)[ l ) +0(e a e 2 ) (23) 

= O ( e a e x ). 

Similarly, exponential decay of and Q a implies (Sbt;. Q) = O (e a ex||CI|fl'i)- The linear term (£q£, Q) is 
zero since C Q C c r a = 0, C Q Q a = -Q ca and = ~Q ca - Lastly, (iV' a (£),Ci) < C1|£||^i by the first 

estimate in Lemma EH Combining the above estimates gives the bound 

ni = o( eo e x |ieii ffl + neii 2 HO- 

By the second inclusion of ©, dcd^Qo. e L 2 (M). Holder’s inequality then implies ||H|| = O (||£||jji). Thus, 
if UWm is sufficiently small, say so that ||-B|| < L, then I + B is invertible and j| (I + B)~ l || < 2. Acting on 
equation m by (I + B) 1 = I — B(I + B) 1 and then llj gives the equation 

( “~ C + V ^ ) =n- a '[X + B(I-B)- 1 X + (I-B)- 1 Y}. 

Using the above estimates of ||B||, [[(/ — H) -1 ||, ||X||, and ||F|| implies 

( a ~ C t V{a) ) = n-'X + O ( C „e x ||^|| H i + Ufm) ■ 

Replacing X by |23|) completes the proof. □ 


5 The Lyapunov Functional 

In the last section we derived dynamical equations for the modulation parameters. These equations contain 
the H 1 ( R) norm of the fluctuation. In this section we begin to prove a bound on £. Recall that the latter 
bound is needed to ensure that u remains close to the manifold of solitary waves M s for long time. 

We employ a Lyapunov argument with Lyapunov function 

T c (i) := A ca (Q ca + £) - A ca {Q ca ) + b'(a){{x - a)Q ca ,g). (24) 

Remark: if f{u) = u 3 , the last term in the Lyapunov functional is not needed; however, apart from compu¬ 
tational complexity, there is no disadvantage in using the above function for this special case as well. 
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Lemma 2. Say u = Q ca + £ is a solution to f7|). where f satisfies {!$ ■ Say e a < 1. If 5'{c) > 0, and 
and IICII// 1 are fess than 1, with ||£||jyi small enough, then 



O ( e a e E + (fiat-xtt + CaCj) ll^llff 1 + Ca^xHCll 


+ HCIIff 1 ) ■ 


(25) 


Proof. Suppressing explicit dependence on x and t, we have by definition 

1 f°° „ 

A ca (u) := Hb(u) — - / u 2 b dx + cP(u). 

" J — OO 

Thus, relations ®, © and © imply that the time derivative of A ca along the solution u is 

+ b" ud x u dx. 


S A “ (u) = 


-cu 2 + b' 


\cu 2 - uf{u) + ^{d x uf + F(u) 


Substituting Q ca + £ for u, manipulating the result using antisymmetry of d x , and collecting appropriate 
terms into b'(a)(C Q ^,d x ((x - a)Q ca )), (AT' a (£), d x [6b(Q ca + £)]), and (A' ca (Q ca ), d x (5b(Q ca + £))) gives the 
relation 

^[A ca(Qca+0 - A ca (Q CQ )] =b'{a)(C Q t,d x ((x - a)Q ca )} +c(Q ca ,f) + {C Q f,d x ( 6 2 bQ ca )) + ||£||| 2 

+ 4 <&'&£) + f - (f'(Qca)£,d x (8bQ) 

+ (Ka(0,d x [5b(Q ca + 0}) + {b"S,d x t) + (A' ca (Q ca ),d x [5b(Q ca +f)}). 


The last term is zero because A' ca (Q ca ) = 0 and since £T Q ca , the quantity c(fi,Q ca ) is also zero. We use 
Lemma 01 assumptions m on the potential, estimates |22l) . and 

\6b'\ < e a e 2 x 

to estimate the size of the time derivative. We also use that Q ca , d x Q ca , d 2 Q ca and f'(Q ca ) are exponentially 
decaying. When e x < 1, higher order terms like (&"£, d x f) are bounded above by lower order terms like {b'f, £). 
Similarly, if ||£|| ff i < 1, then e a £x|l£llffi < This procedure gives the estimate 

j f {A ca (Qca+0 ~ a ca (Q ca )\ =b'(a)(Z,£ Q d x ((x - a)Q ca )) + (N' ca (f), 6bd x 0 

+ O ([c|||£||ffi + ^^IlCllff 1 + eaexllCllffi) • 


Applying the chain rule to the integrand of 



F(Qca+0-HQca)-f(Qca)^ 


2 f(Qca)e 


6b 


dx = 0 


and using the definition of N' ca (f) gives that 


(N' ca (0,Sbd x 0 =(N' ca (0 + -f"(Q c )e,6bd x Q c ) 


j ( F(Q ca + f)- F(Q ca ) - f(Q ca )f - i/'(Qca)£ 2 ) b' dx. 
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The second estimate and the proof of the third estimate of Lemma H of Appendix ^ then imply the bound 
( N L(0> Sba x Q = O (e a e®||Cllffi) ■ Thus, since e a £xM\\ 3 H i < e a e x \\^\\ 2 H1 when ||£|| ff i < 1, we have 


~jj;[A-ca(Qca + £) — A ca (Qca)] — b' (a )(£, Cq8 x ((x — a)Q ca )) + 0 (|c|]]£||#i + e a e^||^||^fi + e a e x ||£||))-i) . (26) 

When f{u) = u 3 , (£,CQd x ((x - a)Q ca )} — 0 since C, ca — cl;c[(x a )Qca]• In this special case the above 
estimate is sufficient for our purposes, but in general, we need to use the corrected Lyapunov functional. 
When £ € C'(R, 7L 1 (R)) fl C 1 (R, H~ 2 ( R)), 6'(a)(£, (.x — a)Q ca ) is continuously differentiable with respect to 
time; 


^ [b'(a){£, (x - a)Q ca )] =d t b'{£, (x - a)Q ca ) + b'(a){£, (x - a)Q ca ) + c6'(a)(£, (a: - a)Q a ) 

+ ab'{a){€, {x - a)Ca) + ab"{a )(£, (x - a)Q ca ), 

where (£,Q C q) = 0 has been used to simplify the derivative. Substituting for using iTHill gives 

^[6'(a)(£, (x - a)Q ca )] = - 6'(a)(£, C Q d x {(x - a)Q ca )) - [a - c + b(a)]b’(a)^\\Q ca \\ 2 L2 + d t b'(£, (x - a)Q ca ) 

+ [a-c + b(a)]b'(a)(d x £, {x - a)Q ca ) + [d - c + 6(a)]6"(a)(£, (x - a)Q ca ) 

+ cb'(a)( £, (x - a)C"a) - &'(a)(£, 8bd x ((x - a)Q ca )) - b'(a)(N' ca (£), d x ((x - a)Q ca )) 
- b'(a)(5 2 bQ ca ,d x ((x - a)Q ca )) + [c - 6(a)]6"(a)(£, (x - a)Q ca ). 

We estimate using the same assumptions used to derive m . If m\w and e x are less than 1, then 

( x - a)Q ca )] = - b'(a){£,C Q d x ((x - a)Q ca )) + O (|a - c + 6(a)|e a e x + |c|e a e x ||£|| H i) 

+ O ( e a4 + ((I + e a)e x + oll^llff 1 + dAi-ll^llffi) ■ 

Adding the above expression to (El gives an upper bound containing |c| and \a — c+ b{a) |. Replacing these 
quantities using the bounds 

|c| = O (e Q e x + £ a c x ||£||z/ 1 + ||£||#i) 

and 

I a — c+ 6(a)| = O (e a e x + IlCHif 1 + ll^llff 1 ) 

from Proposition H and bounding higher order terms by lower order terms gives El- To use the above 
bounds on |c| and \a — c + 6(a)| we must assume ||£||#i is small enough so that Proposition H holds. □ 


6 Spectral Properties of the Hessian 

The Hessian d 2 A ca at Q ca in the L 2 (R) pairing is computed to be the unbounded operator 

C Q :=-d 2 x +c-f'(Q ca ), (27) 

defined on L 2 (R) with domain H 2 { R). We extend this operator to the corresponding complex spaces. 


13 


Proposition 3. The self-adjoint operator Cq has the following properties. 

1- £qCc'„. = 0 and CqQ a = -Q ca . 

2. All eigenvalues of Cq are simple, and Null Cq = Span {£*£}. 

3. Cq has exactly one negative eigenvalue. 

4. The essential spectrum is [c, oo) C R + . 

5. Cq has a finite number of eigenvalues in (—oo,c). 

Proof. Recall that the vectors := ~d x Qca and := d c Q ca are in the Sobolev space H 2 (M.). Thus, 
relations Cq £*£ = 0 and Cq Q a = —Q ca make sense, and are obtained by differentiating A' ca (Q ca ) = 0 with 
respect to a and c. The first relation above proves that Cffj, is a null vector. 

Say (, g € H 2 { R) are linearly independent eigenvectors of Cq with the same eigenvalue. Then, since Cq 
is a second order linear differential operator without a first order derivative, the Wronskian 

W(g,C) = - vdxC 

is a non-zero constant. With g and ( both in H 2 ( R) however, the limit linx,,-^ W(g, £) is zero. This 
contradicts the non vanishing of the Wronskian, and hence all eigenvalues of Co are simple and, in particular, 
Null Cq = Span {<£}. 

Next we prove that the operator Cq has exactly one negative eigenvalue using Sturm-Liouville theory on 
an infinite interval. Recall that the solitary wave Q ca (x) is a differentiable function, symmetric about x = a 
and monotonically decreasing if x > a. This implies that the null vector Cff a , or equivalently, the derivative 
of Qca with respect to x, has exactly one root at x = a. Therefore, by Sturm-Liouville theory, zero is the 
second eigenvalue and there is exactly one negative eigenvalue. 

We use standard methods to compute the essential spectrum. Since the function f(Q C a{x)) is continuous 
and decays to zero at infinity, the bottom of the essential spectrum begins at lim^—xx, (c — f(Q ca 0*0)) = C 
and extends to infinity: a ess (CQ) = [c, oo). Furthermore, the bottom of the essential spectrum is not an 
accumulation point of the discrete spectrum since f(Q C a{x)) decays faster than x~ 2 at infinity. Hence, there 
is at most a finite number of eigenvalues in the interval (—oo, c). For details see 0O1E2- □ 


7 Strict Positivity of the Hessian 

In this section we prove strict positivity of the Hessian Cq on the orthogonal complement to the 2-dimensional 
space df 1 TQ ca M s = Span {Q ca , This result is a crucial ingredient needed to prove the bound on 

the fluctuation £. 

Proposition 4. Assume S'(c) > 0 on I C K.+ . If ^Ti9“ 1 Tq c<i M s , then there is a positive constant p such 
that (CqZ,£) > pU\\ 2 h x. 

Proof. Define X := € JT 1 (R) | £1 ||£|| L 2 = 1}. By the max-min principle, inf X nif»(R)(£ q£, 0 

is attained or is equal to inf cr ess (£Q) = c. If the later holds the proof is complete. In the former case, let g 
be the minimizer. 
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We claim the set of vectors {CtaiCccvV} is an linearly independent set. If they were dependent, then, 
since (, and (fl a are orthogonal, there are non-zero constants a and /3 such that rj = + flCca- Projecting 

this equation onto d and d~ l Q a gives the equations f36'(c) = 0 and aS'(c ) = 0. Thus, the assumption 
<5'(c) > 0 implies 77 = 0. A contradiction since the zero function does not lie in the set X. Note that in 
deriving aS'(c) = 0 we have used that dfl 1 is antisymmetric on the span of Q a since d~ x Q a G A 2 (R). 

By the min-max principle, if 

A 3 := inf {max {(A^,,?} If e V, (IfH^ = 1} | V C H 2 { R), dimP = 3} 

< max{(£ q£, 0I?6 Span {£f, Q a , ;?}} 

is below the essential spectrum, then it is the third eigenvalue counting multiplicity. Let £ = arj + + 7C"q 

where a, j3 and 7 are arbitrary apart from satisfying ||£||/,2 = 1 . Thus, since the third eigenvalue of Cq is 
positive (see Section | 6 j), 

0 < (£ q €,0 = a 2 {C Q r],ri) - 7 V(c) < a 2 (C Q ri,r)) 1 

and hence (£q? 7 , ?y) > 0. The function <r(c) = ( CQr],r 7 ) is continuous since both d~ 1 C*a and d^Cca are 
continuous in A 2 (R) as functions of c. Set g = inf/ <r(c). 

We now improve the result to an IL 1 (R) norm. If we define the constant K(I) := sup 7 ||c — f'(Q C a)\\oc, 
then (£q£,£) > ||d x £|| 2 2 — A"(/)||^|| 2 2 . Adding to this bound the factor A±I 0 f the lower bound (£q£, £) > 
eUWh derived above completes the proof. □ 


8 Bound on the Fluctuation 


We are now ready to prove the bound on the fluctuation. 

Proposition 5. Say e a < 1. Then, for small enough e x < 1 and initial fluctuation ||f(0)||//i < 1, there 
exists a constant C such that the bound 

IlfWIlff 1 = O ( e o + (ea£x ) 2 + e x + 

holds for all times t <T = C (e a e x ) 1 • 


Proof. Lemma |3 implies 


dt 


r dt) 


< c 


( e a e x + ( e a^ x et + e a e 2 x ) ||f || T + e a e x ||£|| T + ||£|h 


for some constant C > 0 where ||£|| T := sup 0<t<T ||£||//i. Integrating over [0, T] gives an upper bound on 
T C (T). A lower bound is obtained by expanding A ca (Q ca + £) to quadratic order then using Proposition 
01 the third estimate of LemmaEland V'(a)(£, ( x — a)Q ca ) = O (ta^Hf Hi/ 1 )- We obtain, after setting all 
non-essential constants to one, 

Ilf IIt ~~ Ilf IIt — e a.e x Ilf || T < r c (T) < |r c (0)| + (e 2 e^ + ( e a e x e t + £a,el) ||f || T + e a^x Ilf Hr + Ilf IIt) ^ 


15 




for all T > 0. Take T = O ^(e a e x ) J ■ Then, under the smallness assumption ||£||#i <C (e 0 e x ) 2 , 

IlCllff 1 = O |T C (0)| 2 + e x + et'j ■ 

The initial value of the Lyapunov functional T c (0) can be bounded by the 7L 1 (R) norm of the initial 
fluctuation ||£(0 )[|jji < Ce 0 (recall that e 0 := infQ oaS M 3 \\uo-QcaWm • Indeed, Taylor expanding A ca (Q ca +£) 
to second order in £ and using the third estimate in Lemma 0 gives |T C (0)| = O (eg + e a e x e 0 ) if e 0 < 1- To 
complete the proof we take e x and eo small enough so that ||£(t)|| i yi is sufficiently small for Lemma|2]to 
hold. □ 

We now prove the main theorem. 

Proof of Theorem 0 By our choice eo < £o, there is a (maximal) time To such that the solution u in 
0 is in U eo for time t < Tq. Hence decomposition 03 with 01, and Proposition 0 are valid for 
the solution u over this time and imply the statements of the main theorem. In particular ||£(f)||#i = 

/ l. i \ ii 

O (e 0 + ( e a e x ) 2 + e x + e t ) for times t < min{T 0 ,T}. Taking eo + (e a e x ) 2 e§ + e x + e t <C e 0 , we must have 

t < T by maximality of the time Tq. □ 


A Estimates of Nonlinear Remainders 


Define 


r°° i 

N ca (0 := - / F(Q C + £) - F(Q c ) - F\Q c )i - -F"(Q C )£ 2 dx 

J — OO " 

and 

Ka(0 ■■= - (f(Qc + 0 - f(Qc) - f\Qc)€) ■ (28) 

Note that N' ca (tf) = d^N ca ((f) under the T 2 (R) pairing. 

Lemma 3. //||£||tfi < 1 and f £ C fc (M) for some k > 3, with f ^ £ T°°(R), then there are positive constant 
C\, C 2 , and C 3 such that 

1- wKaimL^cmh 

2- \\K a (o+y"(Qca)e\\L* <c 2 m 3 H i 

3 . |Ar ca (C)|<C 3 ||a^i- 

Proof. Taylor’s remainder theorem implies 


Ka(0 


k—1 

Y,-J {n HQca)C ~ R{Qca,Z), 
n\ 
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where, since / ^ £ L°°(R), \R(Q C a,f)\ < Cj£| fc . Recall that Q ca is continuous and decays exponentially to 
zero. Together with the assumption that / £ C fc (R), this implies /(")( Q ca ) £ T°°(R) for 2 < n < k — 1. 
Thus, after pulling out the largest constant, 


k 

ll^ca(OIU=<C^||rl|L- 

n=2 

To obtain item0 we use the bound 11C" 11 l 2 < C11 £ 11 , which is obtained from the inequality ||£||l°° < C11 £ 11 #i 
and the assumption that ||£||#i < 1. 

Clearly, slight modification of the above proof gives items 0 and 0 For the latter we use that the 
assumptions on / imply F £ C fc+1 (R) with F'( fc+1 ) £ L°°(R). □ 
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